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Abstract 



The yield of neutrons produced by cosmic ray muons at a shallow depth of 
32 meters of water equivalent has been measured. The Palo Verde neutrino 
detector, containing 11.3 tons of Gd loaded liquid scintillator and 3.5 tons 
of acrylic served as a target. The rate of one and two neutron captures was 
| determined. Modeling the neutron capture efficiency allowed us to deduce the 

total yield of neutrons Y to t = (3.60 ± 0.09 ± 0.31) x 10~ 5 neutrons per muon 
and g/cm 2 . This yield is consistent with previous measurements at similar 
depths. 



I. INTRODUCTION 

In the present work we report the results of a new measurement of the muon-induced 
neutron production at shallow depth, based on the Palo Verde neutrino detector [![]. Even 
though the device was designed for a different purpose, namely the detection of reactor u ei it 
was possible to operate it in a mode suitable for the identification of muon-induced neutrons. 

Neutrons and other hadrons produced by cosmic-ray muons in the earth are an important 
and unavoidable source of background for underground low rate experiments. Knowledge of 
the rates for these processes is an essential part of understanding backgrounds in neutrino 
detectors and other low counting rate experiments. For example, since muon-produced 
fast single or multiple neutrons can mimic the correlated signature of inverse neutron beta 
decay, searches for neutrino oscillations at nuclear reactors must cope with this source of 
background. Other neutrino and proton decay experiments, as well as dark matter searches 
(even though often at greater depth), have to cope with this source of background as well. 
The CDMS experiment, for instance, is searching for cold dark matter 0, and is presently at 
shallow depth; muon-induced neutrons represent a major source of background. Low-energy 
accelerator neutrino oscillation searches are usually performed near the Earth's surface 
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where muon-induced neutrons are a significant source of background. Another example is 
the proposed ORLaND neutrino detector at the Spallation Neutron Source at Oak Ridge 
National Laboratory |J which is also planned to be at shallow depth. 

Despite the importance of the subject, relatively little recent progress has been reported. 
Several measurements of the neutron production rates at various depths conducted in the 
past suggest a smooth dependence of the neutron yield on depth or, equivalently, on the 
average muon energy However, a more recent measurement by the LVD collaboration 

|10| , resulting in a smaller yield, is in disagreement with this simple dependence on depth. 
At shallow depth, where the average muon energy is ~ 10 GeV, the smooth trend was 
confirmed by measurements with relatively small detectors without much shielding against 
the neutrons produced outside them |7|JTT||. In the more recent experiment jlT| , the single 
and double neutron yields were determined separately, and pion production by muons was 
also observed. One drawback of that experiment, possibly present in the other ones as well, is 
the difficulty of distinguishing between the neutrons produced in the detector (the intended 
source of neutrons) and the neutrons produced outside by the muon-induced showers. In 
fact, in Ref. ]TT[ it has been estimated that half of the detected single neutrons originated 
in the hadron showers outside the detector volume. 

The theoretical description of this background process is usually based on the assump- 
tion that the electromagnetic interaction of high-energy muons with matter can be modeled 
by replacing the exchanged virtual photon by "equivalent" real photons, and using known 
photo-nuclear reaction cross sections ||12|| . The analysis is complicated, since in order to re- 
late the theoretical neutron production yield to measurement, the propagation and possible 
cascade multiplication of all reaction products must be understood. (For example, a 7r~ 
produced by a muon will make more neutrons, as does the (n, 2n) reaction, etc.). While the 
smooth variation noted above for the neutron reaction rate versus depth is supported in some 



calculations [13], other approaches come to different conclusions |14|,|15|1. In particular, Ref. 



[IB] was devoted, unlike the others, to the relatively shallow depth relevant for the present 



work. In that work, the calculated neutron yield was smaller than the measured one [|TT 
while the calculated and measured ir + yields agreed with each other. One should also keep in 
mind that while the theoretical description quoted above deals with the muon-nucleus inter- 
action involving the exchange of a virtual photon, neutrons can be also produced by nuclear 
interactions involving real bremsstrahlung photons, and electron/positron pairs created dur- 
ing the passage of muons through matter. (The present experiment cannot separate neutrons 
created by interactions involving virtual photons from those produced by interactions of real 
photons.) 



II. EXPERIMENT 

The measurement of the neutron production rate was performed at 32 ± 3 meter-water- 
equivalent (mwe) depth using the Palo Verde neutrino detector. The detector and its opera- 
tion is described in detail in Ref. |1[. Briefly, the apparatus consists of 66 acrylic cells, each 
900 x 12.7 x 25.4 cm 3 . These cells are filled with liquid scintillator loaded with Gd 0.1% by 
weight, for a total scintillator volume of 11.34 tons. The acrylic material has an aggregate 
mass of 3.48 tons. Muons can thus spall either in the scintillator volume or in the acrylic 
material, and the resulting neutrons cannot be distinguished by source. Neutrons, when 
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moderated to thermal energies, are preferentially captured on Gd, resulting in a ~ 8 MeV 
gamma-ray cascade, the characteristic neutron capture signal. The central active detector 
volume is surrounded by a 1 m thick water shield, and the outermost layer of the detector 
is an active muon veto counter providing 4n coverage. 

Since the apparatus was intended for the detection of reactor neutrinos, it was not 
optimized for the neutron yield measurement. For the present purpose special runs were 
performed without the trigger rejection of the cosmic rays. The throughgoing muon sample 
was selected off-line such that at least two veto hits were recorded. Only muons which 
at the same time went through at least three cells of the central detector were included. 
The delayed neutron capture events, recognized by their energy deposit, are of two kinds: 
The "single bank" events, where only one neutron capture candidate event occurs following 
the veto hit; and the "two bank" events, which have two neutron capture candidate events 
following the muon. (The detector electronics is not capable of recording more than two 
correlated events.) 

The measured quantities are the numbers of single (Ni) and double neutron-capture 
events (N 2 ) associated with muons traversing the central detector. The average path 
length of these muons in the central detector is X (measured in units of g/cm 2 ). Since some 
of the neutrons could have been created in the water shield or in other external detector 
parts, correction factors Q k < 1 are applied to N k , k = 1,2. Finally, the neutron detection 
efficiencies e^/ are introduced. Here, e.g, e^i is the probability that one neutron was created 
and one detected, while ei^ is the probability that two neutrons were created and only one 
detected. If all these quantities were known one could define neutron yields Yi, etc. (per 
muon and g/cm 2 ), where I is the number of neutrons, that are independent of the detector 
properties and obey the relation 

Here e d k aq is the detection livetime correction which in this case depends on the number of 
detected neutron-like events k. 

Obviously, since only two quantities, i\q and N 2 , are measured, only two yields can be 
determined, and the system above cannot be solved without approximations. It is assumed 
further that the contributions from four or more produced neutrons can be neglected, as 
indicated in equation (1). Moreover, as will be shown below, the total neutron yield 

Y tot = Y 1 + 2 • Y 2 + 3 • Y 3 (2) 

is essentially independent of the ratio Y 3 /Y 2 , i.e., on the assumed value of Y 3 , while the 
deduced single and double neutron yields Y\ and Y 2 depend on that ratio significantly. 
Thus, the final results of the present experiment will be expressed as the measurement of 
Ytot- 

Note that the quantities Y k , and thus also Y tot , contain all processes that lead to the 
production of k neutrons by the muon. Thus, if the muon creates a n~ (or any other particle 
except a neutron) which then, in turn, creates k neutrons, all of the neutrons, regardless of 
source, contribute to Y k . This somewhat awkward definition is necessary since the evaluation 
of the efficiencies e k; i is based on a code that tracks just neutron propagation in the detector. 
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The neutron capture events were recorded in two runs, each about half of a day long. 
The raw muon veto rate was about 2 kHz and the rate of muons which went through at 
least three cells of the central detector and caused two or more detectable veto hits in the 
first (second) run was 270(275) Hz. These two runs are essentially equivalent, and the final 
total number of muons was determined to be iV M = 1.42 x 10 7 , very similar in both runs. 

The average path length X of the muons is estimated with a simple ray-tracing Monte 
Carlo simulation, starting from a cos 2 9 zenith angle distribution. The resulting X was 125 
g/cm 2 for the central detector (i.e., the scintillator, the acrylic cells and the small amount of 
other materials (Cu, Fe) in the central detector) and 317 and 62 g/cm 2 for the water shield 
and veto counter, respectively. 

The neutron capture events Nf. were selected using cuts similar to the neutron part of 
the neutrino signal M. For the single neutron events a time cut of 10 - 100 /xs between the 
last muon event and the neutron capture event was used, while for the two neutron events 
only the lower limit of 10 /is was kept, the upper one being irrelevant. The time between 
the two neutron capture events was restricted to 5 - 100 /is. The energy distributions of all 
neutron events had identical shapes within statistics. 

There are several contributions to the livetime correction factors e^ q . The largest one, 
which affects primarily the single neutron events, arises since the initial information is over- 
written when another muon strikes the veto after the single neutron event was recorded and 
before the 450 /js time interval of accepting new events. Combining the correction for this 
effect with other deadtime corrections, the resulting ef aQ = 30.8(27.8)% for the one-neutron 
events in the first(second) runs, and e^ q = 72.4(70.3)% for the two-neutron events. The 
difference in ef a<? between the two runs is caused primarily by the ~ 8% change in the raw 
muon rate. 

The distribution of the time interval between the muon and the single-neutron event is 
shown in the top panel of Fig. 1, while the lower panel shows the same distribution for the 
first of the two neutron events. The fitted curves are of the form 

F(t) = (a ie - t/Tl + a 2 e- t/T2 + a 3 ) • e"'/ 500 ^ . (3) 

Here n = 28.8 ±1.0 [is for the single neutron case while T\ = 11.8 ±0.4 /is (not very far from 
half of the previous Ti, the value one expects for two neutrons) for the double neutron case, 
are the characteristic neutron capture times. Both agree quite well with the Monte Carlo 
(MC) simulation, as seen in Fig. 1. The presence of the time constant r 2 is related to the 
inhomogeneous nature of the detector. Some neutrons enter the acrylic, where there is no 
Gd available for captures, thus prolonging their capture time. Finally, the last term, with 
the fixed time constant of 500 /is, represents the accidental background (500 /xs is the average 
time interval between successive muons). Using such a fit one can subtract the accidental 
background in the relevant time window, and calculate the corresponding correction for the 
described choice of the time cuts. 

In Fig. 2 the distribution of the elapsed time between the two neutron capture events 
is shown, and again fitted to the same functional dependence, Eq. (3). (The distinction 
between the top and bottom panels is explained below.) Here, the time constant T\ = 
18.1 ± 1.0 /is in the top panel is substantially smaller than the T\ obtained by the MC 
simulation {j^ ~ 29 /is, as expected if only two neutrons are involved). This feature 
suggests that the data sample contains a nonnegligible component with more than two 
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neutrons. We will return to this point later. 



III. EFFICIENCIES 

While the evaluation of N^, N\,N2, and X is rather straightforward, the determination 
of the efficiency matrix is somewhat model dependent. This is because, in order to 
calculate e^z, one needs the initial energy and angular distributions of the created neutrons, 
both of which are poorly known. In practice, a number of physically plausible assumptions 
are made about these distributions, and the spread among the resulting efficiencies yields a 
measure of the systematic uncertainty. Once the neutron initial distribution is chosen, the 
Monte Carlo simulation of the neutron transport developed for the neutrino experiment |l ] 
is used. The detector geometry, materials and electromagnetic interactions are simulated 
using GEANT |16[]. Hadronic interactions are simulated by FLUKA |17]], and the low energy 



neutron transport by GCALOR | |18| ]. 

To evaluate the quantities ei 5 i and £2,1 initial single neutrons were distributed randomly 
in position and initial angles through the detector volume. For the initial energy distribution, 
several possibilities were used: the exponential distribution exp~ E ^ 39MeV as proposed by the 
Karmen collaboration |19|], and functions with the power dependence E~ x , 0.5 < x < 2 (see 
|20|). These cover the shape following E~ im suggested by the experience with photo- nuclear 



processes ||21|| . All of these distributions result in similar efficiencies. Using the spread of 
the calculated values as a measure of the systematic error, and taking a simple average as 
the most probable value, one obtains 

ei,i = 0.20 ±0.03, e 2 ,i = 0.007 ±0.004 . (4) 

Note that the probability that a single initial neutron will result in two neutron captures, 
characterized by e 2 ,i, is quite small. 

The determination of two neutron efficiencies is even more difficult, since it depends on 
the energy distribution of both initial neutrons. To evaluate the quantities e^2 and €2,2 two 
neutrons were created at the same random point and with random directions for each of 
them. Two extremes were considered. In the first of them, the total energy of the neutrons 
followed one of the previously described functions used above in the single neutron case. 
This initial energy was distributed randomly between them. (Essentially the same efficiency 
was obtained if the two neutrons shared the initial available energy equally.) The other 
extreme is obtained if both neutrons each have the energy distribution used in the single 
neutron case, i.e., the total neutron energy is on average twice as large. The spread (not 
very large) was again used in the averaging and in the estimate of the systematic error, 
resulting in 

e 1)2 = 0.22 ±0.01, e 2 , 2 = 0.06 ±0.01 . (5) 
Note that there is a sizable probability that only one neutron is detected when two 



were initially created. None of the previous analyses, in particular Ref. |l lfl, took that into 
account. Next, one has to consider the possibility that three neutrons were spalled by the 
muon, but only one or two were detected. Again, using the average and the two extreme 
possibilities to divide the available energy among the three neutrons, one obtains 



5 



€i,3 = 0.19 ±0.01, € 2 ,3 = 0.10 ± 0.01 . (6) 

It is important to realize that these efficiencies are not much smaller than those for initial 
one or two neutrons. Thus, the effect of Y 3 should be considered. (The effect of four and 
more neutron spallation will be neglected, however.) 

The correction factors Qk that exclude neutrons created outside the central detector 
volume (i.e. in the water shield since the effect of the veto is negligible), but captured 
there, must be determined also. To do that, the efficiencies for each passive volume were 
determined using the same Monte Carlo code as for the case of the central detector. Then 
Qk were determined from 

q \X X^Z Yl€k,l) central det. 

kl) central det. 

+ (XEiY l e k , l ) 

water 

resulting in Q 1 = 0.80 ± 0.10 and Q 2 = 0.94 ± 0.07. In order to obtain the above values a 
crude assumption Y\ = 2Y 2 = 2Y 3 , was used. It is important to note that the uncertainties 
of Qk are strongly correlated with the error in the efficiencies. 

The effect of the neutron component of hadronic cascades created outside the veto, and 
coincident with the muon which created them, remains to be determined. The ratio of 
muons making three hits in the veto, thus involving more than one particle, to the prevalent 
case of the two hits was employed as a measure of the frequency of showers. This ratio was 
6.7 ±0.3% and 7.3 ±0.3% in the first and second runs, respectively, i.e., about 7% of muons 
entering the veto were accompanied by a shower. 

To see that the 'triple- veto- hit' events are really different from the standard throughgoing 
muons with just two veto hits, one can form ratios Rk 

]\Tshower 1 ]\Ttotal 

D = iV fc / iV fc /o\ 

k ]\f shower j 'jy total ' V ' 

which represent a quantitative measure of the neutron content of the showers. The measured 
values are R\ = 2.0 ± 0.1 and R 2 = 4.4 ± 0.2, both significantly larger than unity. Thus the 
shower events indeed are richer in neutrons, particularly in the two neutron sample. 

Moreover, the double neutron events with three veto hits have distinctly different time 
structure than the more common ones with just two veto hits. This is shown in the bottom 
part of Fig. 2 for the interevent time, i.e. the elapsed time between the first and second 
neutron captures. The corresponding time constant T\ for the three- veto-hit events is signif- 
icantly smaller than for the two- veto-hit events, showing that the shower events have a large 
multineutron component. A similar effect is present when the capture time of the first neu- 
tron is considered. It is therefore likely that the discrepancies between data and simulation 
in the time dependence of the interevent time interval, noted above, are the consequence of 
a multineutron (k > 2) component in the two bank events. 

Since the shower events contain an unknown number of neutrons created outside the 
detector volume, they are excluded from further consideration. 



IV. RESULTS 

In the first run, the observed numbers of neutron captures, corrected for the random 
background and with the effect of showers subtracted, were Ni = 3916 ± 66, N 2 = 828 ± 29, 
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while in the second run Ni = 3451 ±62, N2 = 829 ±29. The two runs, which were separated 
in time by 10 months, give consistent results when the differences in e daq is taken into 
account, proving that the experiment is stable. 



As a first step, the yields are analysed as in Ref. |11|| , i.e., only the 'diagonal' efficiencies 
6k,k are taken into account, thus 

ysimple = ^ . Qk/{N ^ . x . ^ . e da 9) (g) 

The resulting yields obtained by averaging the two runs are: 

ysimple = , 2M ± . 4( s tat) ± .50(sj/st)) x 10~ 5 neutrons ( 10 ) 

/i ■ g • cm 1 

Y s imp ie = (Q 9g ± o 03 ( sto t) ± O.U(syst)) x 10~ 5 neutrons . 

/i ■ g • cm~ 2 

The exclusion of external showers, characterized by the three- veto-hit events, has resulted 
in the reduction of Yi imple by 14%, and of Y 2 simple by 45%. Thus, presumably due to the 
presence of the sizable shield X ~ 380 g/cm 2 , the effect of the external showers, while still 
clearly present, was reduced compared to the findings of Ref. fLl |, where (in the same units) 
ysimple = 4 3(2.0) x 10- 5 and Y 2 simple = 1.6(0.5) x lO" 5 without (with) the correction for 
external showers. 

However, the proper analysis should include the full efficiency matrix, i.e., the possibility 
that two neutrons were initially produced but only one neutron capture was recorded and 
vice versa, as well as the possibility that three neutrons were originally produced and only 
two or one neutron captures were observed. As pointed out earlier, with only two measured 
capture rates, Ni and iV 2 , it is impossible to deduce all the relevant information without 
further constraints. To avoid these difficulties, and to make the comparison with other 
experiments easier, the total number of neutrons produced per muon was evaluated. This 
quantity, Y tot , defined earlier in Eq. (2), has the further advantage that it is essentially 
independent of the ratio I3/I2 for three to two neutron production, as demonstrated in 
Table I. 

Thus, the final result of the present measurement can be expressed as 

c neutrons 

Y tot = 3.60 ± 0.09 ± 0.31 x 10~ 5 , 11 

/i • g • cm~^ 

where the systematic error is an estimate based on the spread of values in Table I added in 
quadrature to the spread of the evaluated efficiencies. 

All available data on neutron yields are collected and compared in Fig. 3. (It is assumed 
that the other measurements ( fTj-|T0|l) also are really Y tot measurements.) For the measure- 



ments of Ref. [ I I] both results, with and without the correction for external showers, are 
shown. In some of the other measurement the shower contribution was excluded, but it is 
not clear how well this has been done, since at least the results at 25 and 316 mwe were 
obtained with relatively small and unshielded detectors. 

At shallow depth there are now three measurements, with essentially consistent results. 
Clearly, still better and more complete measurements are desirable, in which the full neutron 
multiplicity and energy spectra are determined and the neutrons produced externally are 
reliably identified. 
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TABLES 

TABLE I. Values of the neutron yields Y m , Y U Y 2 , and Y 3 in units of 10" 5 for different 

fi ■ g ■ cm 

assumed ratios YzjYi. 



Assumed T3/T2 


Ytot 


Ti 


Y 2 


T 3 


0. 


3.54 


2.30 


0.62 


0.0 


0.5 


3.60 


2.48 


0.32 


0.16 


1. 


3.62 


2.54 


0.22 


0.22 


2. 


3.64 


2.59 


0.13 


0.26 
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FIGURES 




FIG. 1. Top: The distribution of time between muon hits and the single-neutron capture events, 
fitted to the exponential function, Eq. (3). Bottom: The distribution of time between muon hits 
and the first neutron capture events for the double-neutron events with an exponential fit, Eq. (3). 
In both panels the MC simulation is shown as a histogram, and the fitted time constants t\ and 
T2 are displayed. 
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FIG. 2. Top: The distribution of the elapsed time between the two capture events with the 
two- veto-hit pattern for the double-neutron events with an exponential fit, Eq. (3). Bottom: The 
distribution of the elapsed time between the two capture events for the subset of the double-neutron 
events with three veto hits. In both panels the fitted time constants t\ and T2 are displayed. 
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FIG. 3. Summary of available data on the total neutron yield Y tot versus depth. The points are 
labelled by the corresponding first author. For Ref. [11] two results are shown, with and without 
the effect of external showers. Present work is denoted by an empty circle. 
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